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Abstract 

We introduce a general and accurate method for determining lattice phase shifts and mixing angles, which is applicable 
to arbitrary, non-cubic lattices. Our method combines angular momentum projection, spherical wall boundaries and 
an adjustable auxiliary potential. This allows us to construct radial lattice wave functions and to determine phase shifts 
at arbitrary energies. For coupled partial waves, we use a complex-valued auxiliary potential that breaks time-reversal 
invariance. We benchmark our method using a system of two spin-1/2 particles interacting through a finite-range 
potential with a strong tensor component. We are able to extract phase shifts and mixing angles for all angular 
momenta and energies, with precision greater than that of extant methods. We discuss a wide range of applications 
from nuclear lattice simulations to optical lattice experiments. 
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1. Introduction 

Lattice methods are widely used in studies of quan¬ 
tum few- and many-body problems in nuclear, hadronic, 
and condensed matter systems, see e.g. Refs. CMl. A 
necessary step in such studies is the computation of 
scattering phase shifts and mixing angles from an un¬ 
derlying microscopic lattice Hamiltonian. Remarkably, 
the same problem arises in the context of experiments 
on optical lattices. Several groups have pioneered the 
use of ultracold atoms in optical lattices produced by 
standing laser waves, to emulate the properties of con¬ 
densed matter systems and quantum field theories ||6]- 
[TOl . The basic concept is to tune the interactions of the 
atoms, both with each other and with the optical lattice, 
to reproduce the single-particle properties and particle- 
particle interactions of the “target theory”. Such studies 
often require a more general setup than a simple cubic 
lattice, for instance in the case of the hexagonal Hub¬ 
bard model ini, which closely resembles the physics 
of graphene im and carbon nanotubes ifTSll . Clearly, 
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a robust and accurate method for computing scattering 
parameters on arbitrary lattices is needed. 

For the scattering of particles on a cubic lattice, 
Liischer’s finite-volume method IH uses periodic 
boundary conditions to infer elastic scattering phase 
shifts from energy eigenvalues. The method has been 
widely used in lattice QCD simulations with applica¬ 
tions to different angular momenta GMSl as well as 
partial-wave mixing lIMI . see Ref. El for a recent re¬ 
view. An important advantage of Liischer’s method is 
that periodic boundary conditions are typically already 
used in lattice calculations of nuclear, hadronic, ultra¬ 
cold atomic, and condensed matter systems. Since no 
additional boundary constraints are needed, the method 
is easily applied to a wide class of systems. 

However, Liischer’s method requires that the finite- 
volume energy levels can be accurately determined, 
with errors small compared to the separation between 
adjacent energy levels. This is not practical in cases 
such as nucleus-nucleus scattering, where the separa¬ 
tion between finite-volume energy levels is many orders 
of magnitude smaller than the total energy of the sys¬ 
tem. Fortunately, this problem has been solved using 
an alternative approach called the adiabatic projection 
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method ||2TI - |2^ . There, initial cluster states are evolved 
using Euclidean time projection and used to calculate 
an effective two-cluster Hamiltonian (or transfer ma¬ 
trix). In the limit of large projection time, the spec¬ 
tral properties of the effective two-cluster Hamiltonian 
coincide with those of the original underlying theory. 
This method has been applied to nuclei and ultracold 
atoms, while applications to lattice QCD simulations of 
relativistic hadronic systems are currently being inves¬ 
tigated. 

Since the adiabatic projection method reduces all 
scattering systems to an effective two-cluster lattice 
Hamiltonian, additional boundary conditions can be ap¬ 
plied to the effective lattice Hamiltonian in order to 
compute scattering properties. This opens the door to 
methods more accurate than Liischer’s by removing the 
effects of the periodic boundary conditions, which are 
otherwise a significant source of rotational symmetry 
breaking. One promising approach is to place the parti¬ 
cles in a harmonic oscillator potential and extract phase 
shifts from the energy eigenvalues ESI ESI. Another 
prominent example is the method used in Refs. HSlEll, 
whereby a “spherical wall” is imposed on the relative 
separation between the two scattering particles. Phase 
shifts are then determined using the constraint that the 
wave function vanish at the wall boundary. This method 
has been applied to the two-nucleon problem in lattice 
effective field theory (EFT) lE8l - l3T1l and to lattice simu¬ 
lations of nucleus-nucleus scattering using the adiabatic 
projection method lETIj24ll . 

In spite of such progress in lattice scattering the¬ 
ory, all methods are still lacking in precision, especially 
when partial-wave mixing and high angular momenta 
are concerned. In previous work, numerical approxima¬ 
tions were used for the study of coupled-channel sys¬ 
tems 0. We now describe an extension of the spheri¬ 
cal wall method, which enables an efficient and precise 
determination of two-particle scattering parameters for 
arbitrary energies and angular momenta. We use an¬ 
gular momentum projection and solve the lattice radial 
equation with spherical wall boundaries, supplemented 
by an “auxiliary potential”. We test our method on a 
lattice model with strong tensor interactions that induce 
appreciable partial-wave mixing. We expect our method 
to be applicable in theoretical lattice studies of nuclear, 
hadronic, ultracold atomic, and condensed matter sys¬ 
tems, as well as in the experimental design of optical lat¬ 
tices. While we discuss only non-relativistic wave me¬ 
chanics in our examples here, the extension to relativis¬ 
tic systems simply entails replacing the non-relativistic 
dispersion relation with the relativistic one. 


2. Benchmark system 


We begin with the eigenvalue equation 


-— + V{r,rr,^a,) 


ip = Eip, 


( 1 ) 


where r is the relative displacement, cTi, with i = 1,2, 
are the spins of the two scattering nucleons with rrij^ = 
2p = 938.92 MeV. Following Ref 0, we take 

V = C ^^1 + ^ [3(f • (Ti)(r • (Ta) - cri • (T 2 ] 

^ (^ 275 ) • 

with C = -2.00 MeV and Ro = 2.00 x lO’^ MeV^. 
We only consider states of total intrinsic spin S' = 1. 
The radial equation is 


1 LjL + l) 

2/ir dr'^ 2fir'^ 


Vj{r) 


= Eiljj{r), 


(3) 

where L is the orbital angular momentum and J the total 
angular momentum. The “effective” potential is 


yAr) = c{i + ^-l^ 


exp 



(4) 


for uncoupled channels, and 


Vj{r) = C 



X exp 



2(J-1) 

2J+1 

6y/j(,J+l) 

2J+1 


6y/j{J+l) 

2J+1 

2{J+2) 
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(5) 


for coupled ones. In the continuum, phase shifts and 
mixing angles are obtained by solving Eq. 0 using the 
potentials 0 and 0 with appropriate boundary condi¬ 
tions. 

As rotational symmetry is broken by the lattice, the 
energy eigenstates of Eq. 0 belong to the irreducible 
representations (irreps) Ai, A 2 , E, Ti and T 2 of the 
cubic group SO{3,Z) rather than the full SO(3) ro¬ 
tational group 0|32lE3l. For cubic periodic bound¬ 
ary conditions, as in Liischer’s method iflTll . the cubic 
symmetry remains exact, thus our solutions can still be 
classified by cubic irreps. Nevertheless, the rotational 
symmetry breaking due to the boundaries makes it diffi¬ 
cult to identify states of high angular momentum and to 
extract scattering parameters. In order to remove these 
effects, we impose a hard spherical wall of radius Rw, 


V^V + A9{r-Rw), (6) 
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Figure 1: (Color online) Left panel: Grouping of mesh points accord¬ 
ing to lattice coordinates (p, <p), with lattice spacing a. Right panel: 
Spherical wall radius Rw^ interaction regions I-III as discussed in 
the text and effective potential V/fr) for uncoupled channels with 
^0 = -25 MeV. 


where 6 is the Heaviside step function and A is a (large) 
positive constant, intended to sufficiently suppress the 
wave function beyond Rw (we set A = 10® MeV). 
We take Rw to exceed the range of the interaction, 
such that the boundary is placed in the asymptotic (non¬ 
interacting) region. We also take 2Riy to be less than 
the difference of the box size and the interaction range, 
which ensures that cubic boundary effects remain neg¬ 
ligible. 


3. Angular momentum decomposition 

Let \ f) (g) |S'z) denote a two-body quantum state with 
separation r and 2 -component of total intrinsic spin S^- 
We define radial lattice coordinates {p, ip) by grouping 
equidistant mesh points, as shown in Fig. [T] To con¬ 
struct radial wave functions, we project onto states with 
total angular momentum (J, Jz) in the continuum limit, 
using 

I”) ® I'S'z), (7) 

where the are spherical harmonics with or¬ 
bital angular momentum (L,Lz)- The gg are 

Clebsch-Gordan coefficients. The parentheses around 
J, Jz and L on the left hand side signify that these quan¬ 
tum numbers are not exactly good quantum numbers. 
Note that Eq. 0 is applicable to arbitrary geometries. 
Here, n runs over all lattice points and the “radial shell” 
is given by the integer m. Then, pm is the distance from 
the origin in units of the lattice spacing a, and ()p„,|s| 
picks out all lattice points for which = |fi|. It may be 
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Figure 2: (Color online) Illustration of rotational symmetry breaking 
effects in the Hamiltonian matrix, given in the basis of Eq. jTJ. The 
colors show the magnitude of the matrix elements. To study unphys¬ 
ical mixings, we remove the tensor component of Vj{r). The result¬ 
ing Hamiltonian matrix should ideally be block-diagonal in the S-, 
D- and G-waves etc. Clearly, the matrix elements that cause unphys¬ 
ical mixings are suppressed by several orders of magnitude. In each 
block, the row and column indices represent the radial coordinates of 
the mesh points. For higher partial waves, entire “radial shells” 
vanish due to the angular dependence of the wave function, and such 
redundant rows and columns have been removed. 


practical (especially for non-cubic lattices) to relax this 
condition to include all lattice points with \pm — |ft| | < <5 
for small, positive (5. On the lattice, the \m)g'^’‘ form 
a complete (but non-orthonormal) basis. We therefore 
compute the norm matrix of these states before solving 
for the eigenstates of the lattice Hamiltonian. 

We find that rotational symmetry breaking is almost 
entirely due to the non-zero lattice spacing a. As we 
take a —>^ 0 at fixed Rw, rotational symmetry is ex¬ 
actly restored. The degree of mixing between differ¬ 
ent total angular momenta J and J' is a useful indica¬ 
tor of rotational symmetry breaking. Such effects can 
be interpreted as arising from the non-orthogonality of 
wave functions in different partial waves when their in¬ 
ner product is computed as a sum over discrete lattice 
points. The degree of mixing is difficult to estimate a 
priori, as it depends strongly on the details of the inter¬ 
action. 

Given a simple cubic lattice with a cubic-invariant 
interaction, unphysical J-mixing only occurs between 
cubic irreps of the same type. If the objective is to 
describe a rotationally invariant system on the lattice, 
then we may simply drop all unphysical couplings be¬ 
tween channels with different J. We find that rotational 
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Table 1: Energy levels and differences A (in MeV) with (w/) and without (w/o) unphysical J-mixing matrix elements. In the former case, we 
compute the eigenstates of the lattice Hamiltonian without a spherical harmonic projection. 



Even parity 



Odd parity 



state 

irrep 

w/ 

w/o 

A 

state 

irrep 

w/ 

w/o 

A 


Ti 

0.037 

0.038 

0.001 

VPi 

Ti 

0.917 

0.918 

0.001 

1^112 

E 

2.764 

2.766 

0.002 

VP{F)2 

E 

1.795 

1.796 

0.001 

VD(G)3 

Ti 

3.347 

3.351 

0.004 

VPo 


3.048 

3.053 

0.005 

VGi 

Ai 

6.562 

6.567 

0.005 

VFo 

^2 

4.616 

4.620 

0.004 

VGi 

Ti 

6.624 

6.637 

0.013 


4ll 

4.998 

5.003 

0.005 


symmetry breaking is numerically insignificant at low 
energies in the spherical wall method. Still, it is in¬ 
structive to study the sizes of the unphysical J-mixings. 
For this purpose, we use a simple cubic lattice with 
a = 100 MeV“^ and Rw = 10.02a. In the radial ba¬ 
sis the Hamiltonian matrix becomes nearly block- 
diagonal, with each block corresponding to a specific J. 
The non-block-diagonal elements induce unphysical J- 
mixing. In Table[T] we examine the lowest energy levels 
with and without J-mixing matrix elements. When J- 
mixing is included, we solve directly for the eigenstates 
of the lattice Hamiltonian without a spherical harmonic 
projection. In Fig. we show the Hamiltonian matrix 
elements in the projected basis defined in Eq. (j^. In 
order to focus entirely on unphysical mixings caused 
by rotational symmetry breaking, we have neglected the 
tensor component of Vj{r) in Fig.|^ The magnitude of 
such unphysical mixing matrix elements is found to be 
greatly suppressed. 


4. Auxiliary potential 

We first consider uncoupled channels, where V van¬ 
ishes beyond an “inner” radius R/. A hard wall at Ryi^ 
gives access to discrete energy eigenvalues only, and a 
very large box is needed at low energies. To resolve 
these issues, we define an “outer” radius Rq, between 
Rj and Rw, as shown in Fig. [T] We also introduce a 
Gaussian “auxiliary” potential in region III, 

Vaux(r) = Vo exp [-(r - Rwf/a^] , (8) 

with Rq < r < R\y, where the separation between 
Ro and Rw is chosen such that I4ux is negligible at 
Ro- Note that I4ux vanishes in regions I and II. The 
energy eigenvalues can now be adjusted continuously 
as a function of Vq. In Fig.we show Vj(r) for Vo = 
-25 MeV. 


In order to extract phase shifts, we express '^(r) in 
region II as 

ipir) = Ah']{kr) — Bh'j{kr), (9) 

for Ri < r < Ro, where h'^{kr) and hj{kr) are 
spherical Bessel functions, and k = y/2jlE. The con¬ 
stants A and B can be determined e.g. by a least-squares 
fit in region II. We note that 

B = SA, (10) 

with S = exp(2f(5j), from which 5j can be obtained. 

For coupled channels, ip has two components with 
L = J ± 1. Given Eq. 0, both satisfy the spheri¬ 
cal Bessel equation in region II, and are therefore of 

the form (j^. If we denote A = (Aj_i, and 

B = {Bj-i,Bj+i)'^, the S'-matrix couples channels 
with L = J ± 1. In the Stapp parameterization ll34ll . 


exp(i5j_i) 

exp{i6j+i) 

cos(2ej) 

i sin(2ej) 

i sin(2ej) 

cos(2ej) 

exp(i5j_i) 

exp(/5j+i) 


where ej is the mixing angle. 

When solving S from Eq. ( [T0| as in the uncoupled 
case, we encounter a subtle problem. Eor a simple hard 
wall boundary, only one independent solution per lattice 
energy eigenvalue is obtained. In order to determine 
S unambiguously, two linearly independent vectors A 
and B are needed. In Ref. 0, this problem was cir¬ 
cumvented by taking two eigenfunctions with approxi¬ 
mately the same energy and neglecting their energy dif¬ 
ference. However, such a procedure introduces signifi¬ 
cant uncertainties. 

As the potential (|^ is real and Hermitian, an exact 
time-reversal symmetry results. We now add to Vj{r) 
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Figure 3: (Color online) Phase shifts and mixing angles for J < 4 and 5 = 1. Full, open and “half-open” squai'es con'espond to Vo = 0, 
Vo = —25 MeV and Vq = —20 MeV, respectively. For Vq = —20 MeV, only partial results are shown in order to reduce clutter. Solid lines 
denote continuum results. 


an imaginary component, 


5. Numerical results 


Vj{r) Vj{r) + 


^^aux(r ) 

^^aux(r ) 


( 12 ) 


where C/aux(r) is an arbitrary, real-valued function with 
support in region III only. This leaves Vj (r) Hermitian 
and the energy eigenvalues real, while the time-reversal 
symmetry is broken. Also, ip and ip* are now linearly 
independent and satisfy Eq. Q in regions I and II with 
identical energy eigenvalues. In addition to Eq. ( [T0| , we 
have the conjugate expression. 


A* = SB*, (13) 


and the S'-matrix 

S=[ B A* ][ A B* ]~\ (14) 

from ( [T0] l and ( [T3| l. Phase shifts and mixing angles can 
then be obtained from Eq. Q- Note that the inverse in 
Eq. ( [I4| ) cannot be computed without C/aux(r), since in 
that case A = —B*. We use 


U^n.{r)=Uo6r,ro, ( 15 ) 

for Ro < tq < Rw, where rg is a radial mesh point 
in region III and Uq is an arbitrary real constant. We 
find that the distortion of the energy eigenvalues and ra¬ 
dial wave function introduced by this choice is minimal. 
The same methodology we have applied here for cou¬ 
pled partial waves can also be applied to more general 
problems with different scattering constituents. 


We benchmark our method numerically with the in¬ 
teraction (|^ using a cubic lattice with a = 100 MeV“^ 
(tt/o = 314 MeV), box size 35a, and we take Rj = 
9.02a, Ro = 12.02a, and Rw = 15.02a. Eor all chan¬ 
nels, we use the real auxiliary potential ([^, while for 
coupled channels we add the complex auxiliary poten¬ 
tial ( [T5] | with ?7o = 20.0 MeV and vq ~ Rw- 

In Eig. we show our lattice phase shifts and mix¬ 
ing angles. We compare with continuum results, ob¬ 
tained by solving the Lippmann-Schwinger equation for 
each channel. All our lattice results agree well with the 
continuum ones, from threshold to a relative center-of- 
mass momentum of pcM = k = 140 MeV. We note the 
marked improvement over Ref. lO for the same bench¬ 
mark system. 

6. Application to arbitrary lattices 

While Liischer’s method has been extended to asym¬ 
metric rectangular boxes llTSl . no standard method yet 
exists for an arbitrary lattice. Our method can be used 
to characterize particle-particle interactions on arbitrary 
lattices, in any number of spatial dimensions. This is 
significant for optical lattices, as the lattice geometry 
is then engineered to reproduce the single-particle en¬ 
ergies of a given condensed matter or quantum field 
theoretical system. Anisotropic lattices exhibit more 
breaking of rotational invariance than a simple cubic 
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Figure 4: (Color online) Phase shift for the ^So channel on 

anisotropic rectangular lattices. Circles, triangles, diamonds and 
squares denote results for lattice spacings = 1.2a, = 1.4a, 

az = 3.0a and a^ = 5.0a, respectively. The dashed lines are in¬ 
tended as a guide to the eye. 


lattice does. This is often an essential feature, e.g. in 
the crossover from a three-dimensional system to a lay¬ 
ered two-dimensional one. In Fig. we show the ^5'o 
phase shift on an anisotropic rectangular lattice, where 
the spacing along the z axis, a^, exceeds those along 
the X and y axes, denoted collectively by a. The unit 
cell volume is 100^ MeV“^ in all cases. While we find 
good agreement with the continuum up to — 1.4o, 
this breaks down when becomes comparable to the 
range of the interaction, with increasing deviation at 
high PCM- Such a crossover to two-dimensional behav¬ 
ior can be characterized in terms of mixing between the 
^S'o and (Jz = 0) partial waves, an effect of ro¬ 
tational symmetry breaking. The low-energy particle- 
particle interactions of any lattice system can be simi¬ 
larly described. 


7. Summary and discussion 

We have described a general and systematic method 
for the calculation of scattering parameters on arbitrary 
lattices, which we have benchmarked using a lattice 
model of a finite-range interaction with a strong ten¬ 
sor component. Extensions to more general interactions 
are straightforward. The Coulomb interaction can be 
accounted for by replacing the spherical Bessel func¬ 
tions by Coulomb functions, and by defining the dis¬ 
tance between particles as the minimum distance on a 
periodic lattice. The spherical wall then removes un¬ 
physical boundary effects. When combined with the 


adiabatic projection method, the techniques we have 
discussed can be applied to any scattering system in nu¬ 
clear, hadronic, ultracold atomic, or condensed matter 
physics. We expect our method to be applicable to op¬ 
tical lattice experiments, in addition to its immediate 
usefulness for lattice studies in nuclear, hadronic, and 
condensed matter theory. In fact, the method proposed 
here has already been used to significantly improve the 
adibatic projection method, as detailed in Ref. j^. 
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